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Abstract: We obtain an almost perfect monopole action numerically after abelian 
projection in pure SU(3) lattice QCD. Performing block-spin transformations on the 
dual lattice, the action fixed depends only on a physical scale b. Monopole condensa- 
tion occurs for large b region. The numerical results show that two-point monopole 
interactions are dominant for large b. We next perform the block-spin transformation 
analytically in a simplified case of two-point monopole interactions with a Wilson 
loop on the fine lattice. The perfect operator evaluating the static quark potential 
on the coarse 6-lattice are derived. The monopole partition function can be trans- 
formed into that of the string model. The static potential and the string tension 
are estimated in the string model framework. The rotational invariance of the static 
potential is recovered, but the string tension is a little larger than the physical one. 
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1. Introduction 



The quark confinement is a key problem to understand non-perturbative phenomena 
of QCD. The dual Meissner effect is a promising candidate [1]. 't Hooft suggested 
the idea of abelian projection [2]. SU(N) QCD is reduced to a [/(f)^ 1 abelian 
gauge theory with magnetic monopoles by a partial gauge fixing (abelian projection). 
The dual Meissner effect is caused by the monopole condensation. Monopoles are 
responsible for the confinement as in compact QED [3-8]. 

Monte-Carlo simulations of lattice QCD are a most powerful method to study 
the non-perturbative phenomena. Numerical studies of the abelian projection in the 
maximally abelian (MA) gauge have confirmed the 't Hooft conjecture. In MA gauge 
[9,10], the string tension derived from abelian Wilson loops gives almost the same 
value as that of non-abelian Wilson loop (abelian dominance) [11-13]. Moreover, 
the monopole contribution to the abelian Wilson loops alone reproduces the string 
tension in SU(2) QCD (monopole dominance) [14-16]. The abelian and the monopole 
dominances are seen also in the behavior of the Polyakov loop in T ^ SU(2) and 
SU(3) QCD [17-19]. These results support the 't Hooft conjecture. 

Wilson's idea of a block-spin transformation and a renormalized trajectory (RT) 
are useful when we study the continuum limit on available lattices [20]. The lattice 
action on RT has no lattice artifact and hence it is called quantum perfect action. It 
reproduces the same physical results as in the continuum limit. It is challenging to get 
the perfect lattice action for the infrared region of QCD. For that purpose, we have 
to extract a dynamical variable which plays a dominant role in the infrared region. 
The numerical evidence for the monopole dominance suggests that low-energy QCD 
can be described by an effective action on the dual lattice in terms of a dual quantity 
like monopole currents. It is very interesting to obtain such an effective monopole 
action. Also it is challenging to set the perfect monopole action with the help of block- 
spin transformations for monopole currents. The monopole action can be obtained 
numerically by an inverse Monte-Carlo method [21]. A block-spin transformation on 
the dual lattice can be performed by considering an n-blocked monopole current [22]. 
The detailed study of the effective monopole action in pure SU(2) QCD has been 
done [23-25]. The monopole action determined in SU(2) has the following important 
features: 

1. Two-point interactions are dominant in the infrared region and coupling con- 
stants decrease rapidly as the distance between two monopole currents in- 
creases. 

2. The action fixed seems to satisfy the scaling behavior, that is, it depends only 
on a physical scale b = na((3). This suggests that the action is near to RT. 

3. Monopole condensation seems to occur for large b region from the energy- 
entropy balance. 
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In order to test the validity of the statement that the action is near to RT, we 
need to check the restoration of the continuum rotational invariance. Then, we have 
to determine first the correct forms of physical operators (quantum perfect operators) 
on the blocked lattice. In the above pure SU(2) study, Fujimoto et al. [26] have taken 
the following steps: 

1. Study the renormalization flow on the projected space of two-point monopole 
interactions alone, considering that they are dominant numerically in the in- 
frared region. 

2. The static potential between quark and antiquark can be evaluated by the 
expectation value of the Wilson loop in the continuum limit. Hence the Wilson 
loop can be regarded as the correct operator evaluating the static potential on 
the fine a-lattice also. 

3. Perform the block spin transformation analytically, starting from the two-point 
monopole interactions with the Wilson loop on the fine a-lattice. The quantum 
perfect action and the quantum perfect operator on the coarse 6-lattice can be 
obtained analytically when we take the a — > limit for fixed b = na. 

4. Compare the quantum perfect action composed of two-point monopole interac- 
tions with the effective action numerically determined from the inverse Monte- 
Carlo method. The parameters in the perfect monopole action and the quantum 
perfect operator are fixed then. 

5. The effective action with the quantum perfect operator can be transformed 
into that of the string model. Since the strong-coupling expansion is found 
to work well in the string model, we see that the static potential is estimated 
analytically by the classical part alone and the continuum rotational invariance 
is restored. 

It is not so straightforward to extend these studies to pure SU(3) QCD, since there 
are three monopole currents k^\s) satisfying one constraint Z)a=i &^( s ) = 0- S° 
far, SU(3) monopole action composed of only one kind of monopole current after 
integrating out the other two has been derived in the case of two-point interactions 
[27]. In this case, the same method can be applied as done in pure SU(2) QCD and 
monopole condensation seems to occur from the energy-entropy balance in rather 
strong coupling region. But the scaling seen in the pure SU(2) case was not seen. 
Also it is important to study effective action for two independent monopole currents 
in order to see the characteristic features of pure SU(3) QCD. 

The purpose of this paper is to report the results of the extensive studies of 
SU(3) monopole actions, especially in terms of two independent monopole currents. 
In Section 2 we briefly review the SU(3) monopole current on the lattice and the 
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inverse Monte-Carlo method. The numerical results of the SU(3) monopole action 
are shown in Section 3. In Section 4 we perform a block-spin transformation of 
the monopole current analytically, considering only the infrared dominant two-point 
monopole interactions. In Section 5, transforming the monopole action into that of 
the string model, we calculate the static potential analytically. In Section 6 the string 
tension is estimated from the results of previous Sections. The conclusions are given 
in Section 7. The correspondence between the monopole action and the dual abelian 
Higgs theory (dual Ginzburg-Landau theory) is given in Appendix. 



2. Monopole current and the inverse Monte-Carlo method 

We extract a U{1) 2 link field u^(s) from a SU(3) link field U^(s) after abelian pro- 
jection called MA gauge. In the SU(2) U(l) link field is defined as u^s) = 
diag(e l61 ' 1 ,e l6 ^ ), where 6^(s) = arg[C/ M (s)] aa (a = 1,2). It satisfies det(w M (s)) = 1 
due to Eli 0£ o) = 0. Since ELi ^ ^ in the SU(3) case, the definition of U(l) 2 
link field necessarily becomes more complicated. In this study, we use the definition 
of Ref. [28,29]. 1 The fields transforming as photon fields under U(l) 2 are defined as 
follows: 



0« = arg[f/ M ] aa - -<j>„, M ee ]T arg[f/ M ] 



3' 



e[-7r,7r). (2.1) 

mod27r 



The U(l) 2 link field defined by u M (s) = diag(e l6 ^ \ e ie »\ e 1 ^ ) satisfies det(u M (s)) = I. 

If the U(l) 2 field strength is defined as eJJ = - dj^ (mod27r), then 

£ o e£J = 2nl (I = 0,±1) and is not always zero. When I = +1(-1), is 
redefined. If 9^ is the maximum (minimum) of (O^j, O^J, ©^); ^ ^ s redefined 
as Q^l — 27r(+27r). Others do not change. Then new U(l) 2 field strengths satisfy 

J2 a = 0- Monopole currents are defined as = l/47re^ pCT <9i,Q^ by DeGrand- 
Toussaint(D-T) [7] which satisfies the constraint 

5>« = 0. (2.2) 

a 

We want to get an effective monopole action on the dual lattice integrating out the 
degrees of freedom except for the monopole currents: 

Z = J VU5(X off )A F (U)e- s W (2.3) 
= Jvu[J Vc5(X off )A F (U)e~ s(u) ] (2.4) 
= J Vue~ s ^f (u) (2.5) 



1 Another definition is seen in Ref. [10]. Both definitions are equivalent in the continuum limit. 
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fc(«)ez 

^ fc ( i «) i o d s a fc(°),o e 



we 



S eff (u) 



(2.6) 
(2.7) 



where C/ M (s) = c^(s)m m (s) and X"^ are the off-diagonal part of the following quan- 
tity: 

X(s) = J2[U,(s)A a Ul(s) + Ul(s - jl)K a U^s - p), AJ, (2.8) 
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(2.9) 



A F (U) is the Faddeev-Popov determinant and 5(k^ a \u) is the delta function corre- 
sponding to the D-T definition of the monopole current. 

We can perform the above integration numerically. We create vacuum ensembles 
of monopole currents {k^(s)} using the Monte-Carlo method and the above defi- 
nition of the monopole current. Then we derive the effective monopole action using 
the Swendsen method [21] which is one of the inverse Monte-Carlo methods. 

In order to take the continuum limit, we perform the block-spin transformation 
on the dual lattice by defining an n-blocked monopole current [22]: 



n-1 



#(«)( S W) = £ fcto( ns + (n _ i)£ + % y + jp + ma). 

i,j,m=0 



(2.10) 



If the action numerically obtained satisfies a scaling behavior S[K^ a \n,a(P)] = 
S[K( a \ b = na(f3)\, that is, the action depends only on b, the continuum limit can be 
taken as n — > oo, a(f3) — > for a fixed physical length 6 = na((3). 

The original Swendsen method must be extended due to the conservation law 
of the monopole current d' fl k^(s) = as in the SU(2) case [23]. Let us assume the 
form of the monopole action as S^A;^] = E« GiSi[k^} (Gi is the coupling constant 
for the operator Si[k^]) and define Si[k^] as a part of S^A;^] which contains the 
currents around a specific plaquette (s', /}', v 1 ). When we consider the expectation 
value of some operator Oi[k^}, the following identity holds as in the SU(2) case: 



= (O t [k^]), 

E*w e z <W> o^),o°<[ fc(e) ] ex P(- ^ GA[fc (8) D 
£k(«) e z ^ fc («) 5 s a fc(-),o ex P(- Ej Gj-S'jffcW]) 



EmWgz ^mw,o exp(- Ej GjSj[k'W]) 



(2.11) 
(2.12) 

(2.13) 
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k'l a \s) = fcW(s) + M< a )( W + V+A"W - <W+»"W' - V<WO- ( 2 - 14 ) 

We use this identity to determine the values of the couplings Gi iteratively. Choose 
trial couplings Gi suitably. If Gi are not equal to Gi for all i, then AGi = Gi — Gi 
are estimated from the expansion of Eq. (2.11) up to the first order of AGf 

(Oi - Oi) = (O.S, - 0~S~)AG r (2.15) 

G'i = Gi + AGi are used as the next trial couplings. These procedures continue 
iteratively until the couplings converge. 

3. Numerical results 

Practically, we must truncate the form of the monopole action to derive it numer- 
ically. We know that short-distant and two-point interactions are dominant in the 
SU(2) case. Here, we assume the following form of the SU(3) monopole action: 

S[k^] — Y^Gi(Si[k m ] + Si[k®] + $[fc< 3 >]) (3-1) 

i 

= £ Gi(Si[kM] + Sift 1 *] + Si[-kW - A;( 2 )]). (3.2) 

i 

The Weyl symmetry (the species permutation symmetry of the monopole currents) 
remains after abelian projection. We adopt 27 two-point interactions whose distances 
are up to 3na(f3) [23] and 4, 6-point interactions of the following simple form: 

E{ E (*£°(*)) 2 } 2 , E( E (Wis)) 2 } 3 - (3.3) 

a,s n=— 4 a,s [i=—A 

The lattice size is 48 4 from f3 = 5.6 to f3 — 6.4. After thermalization, 30 configurations 
are used for the average. We determine the lattice spacing a(/3) without using the 
theoretical asymptotic beta function. It is given by the relation a(/3) = \J 'cr(fl) / cr p h 
[30], where a p h is the physical string tension 2 . The results are summarized as follows: 

1. The monopole action for two independent types of monopole currents is ob- 
tained clearly (Figure 1). The qualitative behaviors are the same as in the 
SU(2) monopole action: 

(a) The monopole action has a compact form. The self-interaction 

G\ Eo,v(^"'( s )) 2 is dominant and the coupling constants Gi decrease 
rapidly as the distance between the two monopole currents increases. 

2 If we use <7p{ 2 = 440 [MeV], b = 2a ^ 2 is approximately equal to 0.9 [fm]. 
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Figure 1: The coupling Gi versus b in two-current case. The 4pt couplings are negative 
really. 



(b) The coupling constants have a direction dependence which is expected 
after blocking. Two nearest-neighbor interactions 

G2E a ,s^k^(s)k^(s + fi) and 0^,^^(3)^(3 + 0) 

are quite different for small b region. 

(c) The coupling constants Gi become very small for large b region. 

(d) The simple 4, 6-point interactions become negligibly small for large b > 
2a~^ 2 . Two-point interactions are relatively dominant for large b region. 

(e) The scaling behavior holds well for n — 4, 6, 8 data, if the physical scale 
b = na(/3) is taken in unit of the string tension Ja v h- The action seems 
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Figure 2: The self coupling versus b in the one-current case. 

to be very near to RT on which one can take the continuum limit. 

2. In order to study if monopole condensation occurs by energy-entropy balance, 
we derive the monopole action, considering only one type of the monopole 
current. For simplicity, only two-point interactions are taken into account. The 
scaling was not seen in the previous study where the two loop perturbation 
value was used for a(/3) [27]. When a((3) is fixed by the string tension, the 
scaling is seen beautifully in Figure 2. Since we restrict ourselves to one type 
of the monopole current, the entropy of the monopole current loop is given 
approximately by In 7 x L (L is a length of one long loop). From the previous 
study, we know that only one long loop and some short loops of monopoles 
exist in the vacuum and the value of the action is well approximated as G\ x L. 
Figure 2 shows that the entropy dominates over the energy in the large b region, 
i.e., Gx{b) < In 7. 



4. Perfect action and perfect operator for monopole current 

The monopole action seems to satisfy the scaling behavior. The scaling is the charac- 
teristic of the perfect action. The perfect action reproduces the continuum rotational 
invariance. To test the rotational invariance, we have to determine the form of phys- 
ical operators on the blocked lattice. For that purpose, as done in the SU(2) case, 
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we perform a block-spin transformation from the small a-lattice (a — > 0) to the fi- 
nite b = na lattice analytically, restricting ourselves to a simple case of two-point 
monopole interactions with a monopole Wilson loop. Note that two-point interactions 
are dominant for large b region also in pure SU(3). 

We start from the following action on the a-lattice: 



Z \A = E Vl a) .o^),oexp{- E k<£\s)D (s -s')kP(s') 

k( a )eZ a,s,s' 

+2mJ2N,(s)k£\s)}, (4.1) 



N^s) = E A_1 ( S " s'^wA-W*' + A), d' a S a(3 (s) = Jp(s). (4.2) 

s' 

Here Jp{s) is an electric current around the Wilson loop. Note that in Eq. (4.1), 
the static potential between quark and antiquark with a = 1 is considered. We have 
used the monopole Wilson loop operator as in Ref. [26]. As the surface S a p, we can 
take any open surface with the fixed boundary C, since the monopole Wilson loop 
operator with the closed surface is unity due to 4-dimensional linking number in the 
continuum limit. Here we take the flat surface as S a p. We expand Dq(s) and adopt 
the following first three terms for simplicity: 

D (s - s') = aA-\s - s') + (38 s , s , + 7 A(s - s'). (4.3) 

The above form of the monopole action is derived from the dual Ginzburg-Landau 
(DGL) theory also (for details, see Appendix). 

We perform the block-spin transformation, using the definition of n-blocked 
monopole currents in Eq. (2.10): 

Z[K^,J] = £ ^, 4 Vwoexp{-E kP(s)D (s - s^JV) 

fc( a )gZ a,s,s' 

+2iriY,N„(s)kU(s)} 

s 

n—1 

xll 5 (^£ a) (s (n) ) - E fcJfW + (n- l)jl + iv + jp + ka)) (4.4) 

a i,j,k=0 

= r v 1 r vb r v$ ^ exP {- ^ - 

J - n J - n J ~ n kWez 

+2™ £ N,(s) fcW (s) + i £ 5 (a) WW (*) + < E M«) E *i a) («) 
£ 7 ?V b) )(KJV b) ) - E ^)} (4.5) 

where we have introduced auxiliary fields 0, and 7^ for the constraints J2 a ^ — 
0, & \jk^\s) = and the definition of K^ a \ respectively. Here we have used such nota- 
tions as 2>y = YlaT^l^ ■ There are following identities for integers m^(s), n^(s^) 
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and P M (s): 



exp{2mm {a \s)d^k^\s)} = 1, (4.6) 

n-l 

exp{27rinW( S W)(/fW( a W) - £ = 1, (4.7) 

i,j,k=0 

exp{2mP,(s)Y / ki a \s)} = l. (4.8) 

a 

Then we can change the integral region of 7^, and <\> from the first Brillouin 
zone to the infinite region. Using the Poisson summation formula 

E /[*«] = const, r dF<?> £ exp{2mF^l^}f[F^], (4.9) 

J —OO ,s 



Z[R( a \j] becomes 

POO roo roo roo 

Z[K (a) ,J] = Vj V(j) VB VF J2 

J — OO J OO J — OO J — OO ,/ X „ 

i<- a >ez 

x exp{- £ i* o) («)A>(* - O^'V) + 2™ E^W^W 

a,s,s' s 

+*E(-^ (a) (*) + + 27r/W( S ))FW( a ) 



n-l 



-i E 7i a) (^ (n) ) E ^ a) l- (4-10) 

a,s(") i,J,fe=0 

Writing the lattice spacing explicitly, we get 

n-l 

-ia 4 n E 7i a) (« a s {n) ) E F^(nas {n) + (n - l)a/2 + iaz> + jap + Aw) 

o,s(") i,j,k=0 

= -m 4 E4 a) M^ a) M, (4.11) 



xW(a S )=na 4 E7i a) (na S W) 



„(n) 



n— 1 

x5(nas^ + (n — l)a — as^) 

II E 5 ( nas S n) + /a - a ^)- ( 4 - 12 ) 

i(^) J=0 



Integrating out and S^ a \ we find 

/OO _ _ /"CO 

T? 7 exp{io 4 E li a) (bs in) )Kl a) (bs {n) )} / T?0 
a, s («) ' J ~°° 

x exp{-ia 8 E ^ (a) M V(as - as')V^(as')} 



s,s' 



x E exp{-n 2 a 8 l ^( as ) A A as - as ') l i a \ as> ) 

al( a )eZ a,s,s' 

+7r 2 a 8 E V^ a \as)A llu {as - as')l^ (as)} , (4.13) 

a,s,s' 
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where 

Vf\as) = Xj?>(as) - ^(as) - 2nN li (as)6 a , 1 , (4.14) 
A^as - as') = j 5, u - ^^D^\as - as'). (4.15) 

Here we consider = alone, since we finally take the a — > limit in the sum with 
respect to al^ E Z. We change the variables from 7^) to 7'^: 

7 '(!) = 7 (2) _ 7 (3) ? y(2) = 7 (3) _ 7 (1) ; y(3) = 7 (1) _ 7 (2)_ 

Integrating out 7 /<a ' ) and M , we get the monopole action in terms of as follows: 

2 

Z[J] = exp{ — it 2 a 8 N lx (as)A IJil/ (as — as')N u (as') 

+ ^ 2 o 8 ]T S M (teW)^(teW-^)S v (^ B ))}Z^[J], (4.17) 

s (n) s /(n) 



where 



1 71—1 



^(fe (n) ) = lim —a 8 ^ 5(nas^ n) + (n - l)a - as M ) JJ ^ 5(nasr + Ia- a Si ) 

n^oo ^ S,S' ,V i^^) ^=0 

x I V - ^^7} £>o V - as')^K), (4-18) 

A >^ bs (n)_ bs >(n) ) = 

n— 1 

lim —a 8 X! <K nas ^ n ' ) + (^ — l)a — as^) J] ^ <5(nasj n ' ) + la — as*) 



x ^ 5™ - 



-j/^as-as') 



n-l 

x5(na4 n) + (n- 1) a — as,, 

n E^ as t )+ ja - as P- ( 4 - 19 ) 

ZmmjJ] is the dynamical monopole part: 

xexp{-6 8 £ ^"H^M^^ -bs' in) )K^\bs' {n) ) 

0jS (n) ja '(n) 

+2«6 8 £ S A1 (6s (n) )A^ 1 (6s (n) -^^^(fts'W)}. (4.20) 

s (n) jS /(n) 

The interactions of the perfect action on the 6-lattice in Eq. (4.20) depend on direc- 
tions. This is consistent with the numerical data. 

The spectrum of K^(bs^) is found to be equivalent to that in the continuum 
theory as discussed in the SU(2) case [26]. 
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5. String representation and rotational invariance 

When we transform Z mon [J] in Eq. (4.20) into the string representation as in SU(2) 
[31,32], we can estimate the static potential analytically. 
First we change the variables as follows: 

K 0) =i W_ j (S) ) KW=jM-jW, K®=jM-jW. (5.1) 

The conservation laws of monopole currents & = leads us to d'^j^ = d'^j^ = 
fffjjp- The conditions are expressed as 

f 2?^(i;v (o) )exp{<E4 a) (*)^ (a) («)}- (5-2) 

J n a a,s 

Then Z mon [J] is reduced to the following: 
Z mon [J]= r T><p6(E<pV) E w{-2E#Hs)A'r(*-^A a) W 

7r a j(")eZ a,s,s' 

+2 E ti\s)A;l{s-s>)^{s)+iY.ti\s)d^ {a) {s) 

a<b,s,s' a,s 

+2mY / B,(s)A'; 1 u (s - s>)(tf\ S ') -^V))} (5.3) 

s,s' 

/oo v fir ^ 1 v . 

DC / zv*(Ev (0) ) E exp{--E^,(^V(^-^)C,( S ') 

+2«E^( S )A- 1 ( S - S ')0?V) - # (*'))}, (5-4) 

s,s' 

where auxiliary fields = (C^, C®) are introduced and e* a are the SU(3) root 
vectors: ei = (1,0), e* 2 = (-1/2, -y/Z/2), e 3 = (-1/2, >/3/2). Since Ea^ (a) = and 
e a = 0, we get an identity 

ex P{"^ E M«) E(W a) («) + ^ • = 1- (5-5) 

The Poisson summation formula 

E exp{27riE W } + %EW = E II^ o) + t-i a) ) (5-6) 

,(«) ez a a j( a) eZ a 

gives us 

roo 



/oo _ 
^ E ^),oexp{2«E^ a) (^)^ ) (^)}/(^ a) ) = 

const. /_~D0 E /O?-^)- (5-7) 
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Then Z mon [J] becomes 

/OO POO v PIT 

VF VC V^5(Zv {a) ) E *e.i«,o 

J -°° a l{a)eZ 

x exp{-^ E ct^A'^s - s')^ u (s') 

s,s' 

+*E^ a) (*)(^ (a) (*) + e a ■ C^s) + 2^(8)) 

a,s 

+2ttE(^ 2) («) - if '(a))^* - s'W*')}- (5.8) 

s' 

Here we perform the Berezinski-Kosterlitz-Thouless (BKT) transformation [33]: 

<9 M y^( S ) + 27rf«(a) = d^}{s) - 2ttE^A- 1 ( S - *>&V), (5-9) 

s' 

where <pW (s) = y? (a) (s) - 2tt £ s , A" 1 (s - s')d' v s^ (s') + 2vrr( a ) (s) is non-compact. The 
plaquette variable aj%) satisfies a conservation law d[ a a^l(s) = and a constraint 
Y.a^} — due to Y.a^jf 1 — 0- Using the condition d'^B^ = 0, we integrate out Fjf-\ 
and (p$. We get the string representation: 

Z str [J] = eM-^ 2 T, B ^) A> ,l(s-s')B u (s')} £ 8 ^),o 



3 

2tt_ 2 
3 



+-T- E(<$(«) " ^RA^s - (5.10) 



s.s' 



- s ') = E - Si)A~ 2 ( Sl - s')- (5.H) 

Sl 

Since the numerical results show that the couplings of the monopole action are 
weak for large b region, we can use the strong coupling expansion in the dual string 
model. Then we see that the quantum fluctuations can be neglected as in SU(2) case 
[26]. The classical part in Eq. (5.10) cancels the second classical term in Eq. (4.17). 
As a result, the classical part of the expectation value of the Wilson loop is reduced 
to 

2 r°° 

(W(C)) cl = exp{--7r 2 / d 4 xd 4 y%(x) A^x - y)N u (y)}. (5.12) 

Since the classical part is written in the continuum form, the continuum rotational 
invariance is trivial. The static potential V(Ib, 0, 0) and V(Ib, lb, 0) can be evaluated 
as in SU(2). We take the following plaquette sources as S a p(x) for V(Ib, 0, 0): 

S a <p(x) = S*iS^0(xiW(Ib - Xl )9(x 4 )9(Tb - x 4 )5(x 2 )5(x 3 ) (5.13) 
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and for V(Ib, lb, 0): 



S a< ^(x) = (5 al 5p 4 + 8 a 28fu)8(x z )6(x4)6(Tb - x 4 ) 

x6( Xl )6(Ib - x 1 )9(x 2 )e(Ib - x 2 )5(x 1 - x 2 ), (5.14) 

respectively. Using the following formula 



2 



l im f^LX = 7rT5 ( a ) ? (5.15) 



we get the static potentials 



V(lb,0,0) = - lim ±-]n{W(Ib,0,0,Tb)) d 

1 ^oo ± 



> 2vr 2 /r ^ 2 f d 2 p 



i^oo 3 v 'J 2tt 2 



AD 



(0,p 2 ,P3,0) 



= — ln w (5 - 16) 

V(Ib,Ib,0) = ^™ Ib ln ^ (5.17) 
3 m 2 

where n{m\ —m\) = 7 -1 , + 7713 = /?/7, m\m\ = a/7. The static potential has the 
linear term alone and the rotational invariance is recovered completely. The string 
tension is evaluated as 

™, m 1 

<Td = — m — . (5.18) 
3 m 2 

This is consistent with the results [34]. The mj" 1 and the m^ 1 could be regarded as 
the coherence and the penetration lengths in Type-2 superconductor. 



6. Estimate of the string tension 

Ths string tension in Eq. (5.18) is written by the parameters of the monopole action 
taken on the a-lattice. The parameters can be determined by comparing the theoret- 
ical perfect action in Eq. (4.20) with the numerical results of the monopole action on 
the 6-lattice. A'~l in Eq. (4.19) fixed the gauge is calculated in ref. [26]. The results 
of the comparison are seen in Figure 3. The agreement is rather good. The value of 
the string tension is \Ja c i/ffph ~ 1-22 for large b as in Figure 4. The result is not so 
bad, although the discrepancy is not negligible. 

In order to study whether the difference is due to the ambiguity of the fit of 
Do, we try to estimate the string tension without determining the parameters on the 
a-lattice. The following monopole Wilson loop on the 6-lattice is considered: 

WC)) = i £ ^^/^^expi-E^H^^^^-^)^^') 

K(°-)€Z M M a,s,s' 

+2mY / N,(s)Kl(s)}, (6.1) 
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Figure 3: The coupling d versus distance. The data from (i) perfect action (diamonds) 
and (ii) Swendsen method (circles). The values of dimensionless k, mi and m,2 are shown. 



where the naive monopole Wilson loop operator (A. 15) on the coarse b- lattice is used 
and D~l corresponds to A'~. Transforming Eq. (6.1) to the string representation and 
neglecting the quantum fluctuations as in the previous section, we get the classical 
part of the Wilson loop as follows: 

(W(C)) d = exp{-^ 2 J2 N ^ s ) D ^( s - sW* 7 )}- (6-2) 

6 s,s> 

It is found theoretically that the static potential V(Ib, 0, 0) evaluated by Eq. (6.2) 
agrees with Eq. (5.16) in the / — > oo case. Using Eq. (6.2), the string tension becomes 

7r 2 d 2 k 1 
ad = ~6 J-n (2tt) 2 (sin 2 f + sin 2 f ) 2 
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Figure 4: The string tension versus b. The data from (i) Eq. (5.18) (diamonds) and (ii) 
Eq. (6.3) (circles). 



x [sin 2 ^D(0, k 2 , k 3 , 0; 2) + sin 2 ^D(0, k 2 , fc 3 , 0; 3)], 



D^{s - s') = 5, v / -——D(ki, k 2 , h, k 4 ; jj)e lk ^ s \ 



(6.3) 



(6.4) 



» v (2tt) 4 " 

We estimate the string tension from Eq. (6.3) with the numerical data of the monopole 
action on the 6-lattice. 

The results are seen in Figure 4. The results are not so different from those 
of the first method for large b. The 6-independence holds roughly for b > 2 and 



Gd/vph ~ 1-27 for large b. The difference from the physical value is also seen in 
this method. 

Using the same method as in SU(2) [26], the effect of the quantum fluctuation 
can be evaluated. The order of the correction for the string tension becomes 



--exp{-5n(0)6 2 }, 



(6.5) 



where 11(0) represents the coupling of the self-interaction in the string model. For 
example, the value of Eq. (6.5) is —2.6 x 10~ 10 for b = 3.27. The correction is very 
small in large b region. 

In order to clarify the origin of the difference, we have to study the effects of 
systematic errors carefully. 
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As a systematic error, the difference between D-T monopoles used in the sim- 
ulation and the monopoles in the analytical calculations may be important. The 
magnetic charge of the latter ranges from — oo to oo, while D-T monopoles take only 
the restricted values of magnetic charges. This is under investigation. 

The glueball spectrum can be evaluated similarly as in SU(2) QCD [26]. The 
lightest ++ glueball mass is determined as M ++ = 2m 2 . For b = 4.23, 2m 2 / v /o^ = 
3.90. The result is not so different from 3.64 ± 0.09 [35]. 

7. Conclusions 

We have studied an effective monopole action in pure SU(3) lattice QCD. The fol- 
lowing results are obtained: 

1. We have found that the SU(3) monopole action can be derived clearly for 
various n-blocked monopoles in MA gauge numerically, using the extended 
Swendsen method in the two-current and the one-current cases. We perform 
n = 2 ~ 8 step block-spin transformations on the 48 4 lattice. In the two-current 
case, we obtain an almost perfect lattice monopole action for the infrared region 
of SU(3) QCD, since the action seems to depend on a physical scale b alone. 
The simple 4, 6-point interactions become negligibly small as compared with 
two-point interactions for large b region as in SU(2). Thus two-point monopole 
interactions are dominant in the infrared region. In the one-current case, if the 
physical scale b is taken in unit of the string tension y/cr p h, the scaling which was 
not seen in the previous study [27] holds beautifully. Monopole condensation 
occurs in the large b region from the energy-entropy balance. 

2. The perfect action satisfies the continuum rotational invariance. When we try 
to study the restoration of the rotational invariance of the static potential, 
we have to know the correct form of the perfect operator. For that purpose, 
a block-spin transformation from the small a-lattice (a — > 0) to the finite 
b = na{P) lattice is performed analytically in a simple restricted case of a 
quadratic monopole action as in SU(2) case. The perfect monopole action and 
the perfect operator evaluating the static potential between electric charges on 
the 6-lattice are derived. The quadratic SU(3) monopole action is simple, but 
it corresponds to the London limit of the DGL theory which is non-trivially 
interacting field theory. 

3. The SU(3) monopole action can be transformed into the string model. There 
are three strings satisfying one constraint J2a a ^ — 0- This is consistent 
with the results in the continuum limit [36,37]. Since the monopole interactions 
are weak for large b region, the dual string interactions are strong. Using the 
strong coupling expansion, the static potential in the long distance is calculated 
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by the classical part alone analytically. The static potential has a linear term 
alone and the restoration of the rotational invariance can be seen explicitly. The 
string tension is estimated from the numerical data of the monopole action. The 
result is rather good, but the difference from the physical string tension is not 
negligible. The same thing happens in SU(2) case [26]. It seems that there still 
exist some systematic errors. 
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A. Appendix 

In SU(2), it is already known that the monopole action can be derived from the 
dual abelian higgs model on the fine a-lattice [8,32,38]. We show how to derive the 
monopole action from the DGL theory on the lattice in SU(3). In DGL, = 
(AJ,AJ) are used as U(lf photon fields, where = A^X/2 (X 1 are Gell-Mann 
matrices). The magnetic charges rft are distributed on the SU(3) root lattice: rft = 
9mY^a=iCa^a, (,a £ Z , g m = 4n / g (g is the SU(3) coupling constant). The DGL 
Lagrangian in the continuum limit [34] is 

Ugl = \(dfip - drf,) 2 + ]T(|(^ + ig m e a ■ SM a) \ 2 + K\<P {a) \ 2 - v 2 ) 2 ), (A.I) 

4 a=l 

where ~5 ^ = {C Z ,CT) are dual U{1) 2 gauge fields and 0^ = exp(i</?^) are 
monopole fields with a constraint J2 a V 9 ^ = 0- 

In the below, we use the notations of differential forms on the lattice [39]. The 
lattice DGL action [40] becomes 

s DGL [^A {a) ] = ^ll^ll 2 -7E(^ a) *^ a) W2 A + M 

Um a,s,n 

+AE(^ a) vf ) - l ? + E^ a) V<°\ (A.2) 

a,s a,s 

where = exp(ie a ■ ). We modify the partition function of the DGL model using 
the Villain approximation e QCOS ^ — > J2iez e~f^ +27r ^ 2 . In DGL, the summations of 
appear due to 2tt periodicity of cos(d(p ( - a ^ + e a ■ C). Then there is a constraint 
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J2 a ^ = 0) since J^aidp^ + £ a ' C) — 0- The partition function of the DGL theory 
is given by 

DC / TV5(Ev? (a) ) / Z>P" 2 E *e.iO>,o 
x exp{--U|^|| 2 - 7 E Pi a V2 A (^ (a) + 6a • £ + 27rf<«% 

-ae^W-^-E^) 2 }- ( A - 3 ) 



Inserting 

i = { n (47Pi a VS,)- 1/2 | /°° ^^ex P{ - e -4^r 

x(F^(s) - 2i 7 pWpg A (V a) + S ' ^ + 27r/W) s ,,) 2 } (A.4) 

and using the relation in Eq. (5.7), we integrate out 0, and The partition 
function is written by the variables in the r.h.s. of Eq. (5.7) and p( a \ From the 
integrals of (p( a \ the relations = 8j^> = Sj® appear. Here monopole currents 
are introduced as follows: 

= jW-j®, = = - j (2) . (A.5) 

They satisfy the conservation law ^/c*"- 1 = and the constraint J2 a k^ = 0. The 
partition function can be rewritten as follows: 

Z mon [k^}= E ^(«),o^ afcM ,oexp{-5l n [A;( a )]-^ on [A;( a )]}, (A.6) 



SZon[k {a) } = 9 -fj:(k^, A" 1 ^)), (A.7) 

^ a=l 



l ^E a pj a Vg A (^(^))^ 



S^k^] = -M / Dp- 2 exp{--E 



^ ' s,n 2^a<bPs Ps+p, 

-AE((pi a) ) 2 -l) 2 "E(pi a) ) 2 }}- (A.8) 

a,s a,s 

The integration of p^ in Eq. (A.8) can't be performed exactly. In the London limit 
(A — > oo), p( a ) are fixed to unity. S^ lon [k^] becomes the quadratic self interaction: 

^[^l^^EH^II 2 - (A.9) 
The monopole action in Eq. (4.1) 

Smon[k {a) ] =J2(k {a \D k^), D = aA- 1 +(3 + 1 A (A.10) 
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corresponds to the modified London limit of DGL: 

X {a) = d(p^ +e a -~d + 2nl( a \ (A.ll) 

The Wilson loop is also calculated similarly. For simplicity, we consider the case 
of the London limit. Since the DGL theory is the dual theory for a color electric 
U(l) 2 charge, we consider the 't Hooft loop: 

"j POO , P7T _ _ 

M = ^- vc / xv5(I> (a) ) E ^k-),o 

ZDGL J-oo J-n a l(a)&z 

x exp{--^| \dd - A^;S\\ 2 - 7 E W a) + ■ ~5 + || 2 }, (A.12) 

Um a 

where S is a surface with the fixed boundary C: 8S — j, where j is the unit current 
on the loop C. Q c are color electric U{\) 2 charges of quarks with colors c (c=R,G,B): 

3 « =(«•«)={& 275) • (4 ^)-( --^)}- « A13 » 

For c=R, we get 

(Hr(C)) = -^- E ^),o^ afe ( ),oexp{-^E(^ (a) > A " lfc(a) ) 

! Y,\\k (a) \\ 2 + ^CS, A-'dk^-^Cj, A" 1 *;)}- (A.14) 



127 V" " " ~"" v ~" ' 

The monopole Wilson loop operators are 

e 2ni(*S, A-Mfef 3 )) e 2ni(*S, A^dk^) g27ri(* S, A" 1 ^ 1 )) ^ 

for c=R,G,B, respectively. When *S is a closed surface, (*£, A _1 <i£; (a) ) in Eq. (A. 15) 
is the 4-dimensional linking number between the closed surface *S and the closed loop 
k^ a \ Hence the monopole Wilson loop is independent of the choice of the surface S. 
This is very important. 

Let us consider the string representation of Eq. (A.12). Introducing the following 
U(l) 2 electric charges: 

s c = ( a « sf\ S W) = {(1, -1, 0), (-1, 0, 1), (0, 1,-1)} (A.16) 

for c=R,G,B respectively [36], the string representation can be derived as in Section 
5: 



(H C (C)) = -L E ^ ^p{-47r 2 7 {E(^ a) , (A + m 2 )-\ 



+AE0?>, (A + m 2 )- 1 ^)}}, (A.17) 
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where = s^j are electric U(l) 2 currents and m 2 = Sjg^ is a mass of a dual 
gauge field in DGL. The two-forms = si a) S + a (a) satisfy 5*a^ = with the 
constraint J2 a a ^ — 0- Thus we obtain the hadronic string model in the SU(3) lattice 
QCD. 
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